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ABSTRACT 

The universe should be dark at energies exceeding ~ 5 x 10 19 eV. This simple 
but solid prediction of our best known particle physics is not confirmed by 
observations, that seem to suggest a quite different picture. Numerous events 
have in fact been detected in this energy region, with spectra and anisotropy 
features that defy many conventional and unconventional explanations. Is there 
a problem with known physics or is this a result of astrophysical uncertainties? 
Here we try to answer these questions, in the light of present observations, while 
discussing which information future observations may provide on this puzzling 
issue. 



1 Introduction 

One of the major goals of cosmic ray physics has always been the discovery and 
understanding of the end of the cosmic ray spectrum. Until the end of the '60s, 



this search was mainly aimed to understand the limits to the acceleration pro- 
cesses and nature of the sources responsible for the production of the particles 
with the highest energies. However, after the discovery of the cosmic microwave 
background (CMB), it became soon clear that the observed spectrum of the 
cosmic radiation had to be cut off at a "natural" energy, even if an ideal class 
of sources existed, able to accelerate particles to infinite energy. In fact, if the 
sources are distributed homogenously in the universe, the photopion produc- 
tion in the scattering of particles off the CMB photons imply a cutoff in the 
observed spectrum of the cosmic rays at an energy Eqzk ~ (4 — 5) x 10 19 eV, 
close to the kinematic threshold for that process 0) . This cutoff has become 
known as the GZK cutoff and particles with E > Eqzk are usually named 
ultra- high energy cosmic rays (UHECRs). 

Several experiments have been operating to detect the flux of UHECRs, 
starting with Volcano Ranch H) and continuing with Haverah Park H) and 
Yakutsk @) to the more recent experiments like AGASA , Fly's Eye 

|, and HiReg |l2|) Thc search for thc GZK cuto ff ? instead of confirm- 
ing the simple picture illustrated above, has provided stronger and stronger 
evidence for the existence of events corresponding to energies well in excess of 
Egzk- the GZK cutoff has not been found. 

This problem hides many issues on plasma physics, particle physics and 
astrophysics, that in their whole represent the puzzle of UHECRs. In section 
|^ thc present status of observations of UHECRs is summarized; in section |^ 
some caveats in the arguments that are often used to address the problem of 
UHECRs are considered, for the purpose of stating the problem in a clear way. 



In sections 4.1 and 4.2 the bottom-up and top-down models of UHECR origin 



are summarized. In section |] some speculations are discussed of new physics 
scenarios that might play a role not only for the explanation of UHECRs, but 
also for the understanding of other current puzzles in high energy astrophysics. 
Conclusions are reported in section ^|. 



2 Observations 



The cosmic ray spectrum is measured from fractions of GeV to a (current) 
maximum energy of 3 x 10 20 eV. The spectrum above a few GeV and up to 
~ 10 15 eV (the knee) is measured to be a power law with slope ~ 2.7, while 
at higher energies and up to ~ I0 19 eV (the ankle) the spectrum has a steeper 



slope, of ~ 3.1. At energy larger than 10 19 eV a flattening seems to be present. 

The statistics of events is changing continuously: the latest analysis of the 
"all experiments" statistics was carried out in lil) where 92 events were found 
above 4 x 10 19 eV. 47 events were detected by the AGASA experiment. A more 
recent analysis H) of the AGASA data, carried out expanding the acceptance 
angle to ~ 60°, has increased the number of events in this energy region to 59. 

In l> the directions of arrival of the AGASA events (with zenith angle 
smaller than 45°) above 4 x 10 19 eV were studied in detail: no appreciable 
departure from isotropy was found, with the exception of a few small scale 
anisotropies in the form of doublets and triplets of events within an angular 
scale comparable with the angular resolution of the experiments (~ 2.5° for 
AGASA). This analysis was repeated in llH) for the whole sample of events 
above 4 x 10 19 eV, and a total of 12 doublets and 3 triplets were found within 
~ 3° angular scales. The attempt to associate these multiplets with different 
types of local astrophysical sources possibly clustered in the local supercluster 
did not provide evidence in that direction 0). 

Recently, the AGASA collaboration reported on the study of the small 
scale anisotropies in the extended sample of events with zenith angle < 60°: 5 
doublets (chance probability ~ 0.1%) and 1 triplet (chance probability ~ 1%) 
were found. 

The information available on the composition of cosmic rays at the highest 
energies is quite poor. A study of the shower development was possible only 
for the Fly's Eye event 0) and disfavors a photon primary \i3/ . A reliable 
analysis of the composition is however possible only on statistical basis, because 
of the large fluctuations in the shower development at fixed type of primary 
particle. 

The Fly's Eye collaboration reports of a predominantly heavy composition 
at 3 x 10 17 eV, with a smooth transition to light composition at ~ 10 19 eV. 
This trend was later not confirmed by AGASA H' HI) . 

Recently in ref. the data of the Haverah Park experiment on highly 
inclined events were re-analyzed: this new analysis results in no more than 30% 
of the events with energy above 10 19 eV being consistent with photons or iron 
(at 95% confidence level) and no more than 55% of events being photons above 
4 x 10 19 eV. 

Recently a new mass of data has been presented by the HiRes experiment 



\1B) . Only two events with energy above 10 eV have been detected by this 
experiment insofar, compatible with the presence of a GZK cutoff. This dis- 
crepancy with the results of several years of AGASA operation needs further 
investigation. Several sistematics have been identified that might considerably 
affect the determination of the energies and fluxes of fluorescence experiments 
versus the ground array techniques H). These issues will not be discussed 
further in the present paper. 

3 The GZK cutoff: how serious is its absence? 

The puzzle of UHECRs can be summarized in the following points: 

• The production problem: the generation of particles of energy > 10 20 eV 
requires an excellent accelerator, or some new piece of physics that allows 
the production of these particles in a non-acceleration scenario. 

• The large scale isotropy: observations show a remarkable large scale isotropy 
of the arrival directions of UHECRs, with no correlation with local struc- 
tures (e.g. galactic disk, local supercluster, local group). 

• The small scale anisotropy: the small (degree) scale anisotropics, if con- 
firmed by further upcoming experiments, would represent an extremely 
strong constraint on the type of sources of UHECRs and on magnetic 
fields in the propagation volume. 

• The GZK feature: the GZK cutoff is mainly a geometrical effect: the 
number of sources within a distance that equals the pathlength for pho- 
topion production is far less than the sources that contribute lower energy 
particles, having much larger pathlength (comparable with the size of the 
universe). The crucial point is that the cutoff is present even if plausible 
nearby UHECR engines are identified. 

• The composition: it is crucial to determine the type of particles that 
generate the events at ultra-high energies. The composition can be re- 
ally considered a smoking gun either in favor or against whole classes of 
models. 

The five points listed above are most likely an oversimplification of the 
problem: some other issues could be added to the list, such as the spectrum, 



but at least at present this cannot be considered as a severe constraint. On the 
other hand, specific models make specific predictions on the spectral shape, so 
that when the results of future observations will be available, this information 
will allow a strong discrimination among different explanations for the origin of 
UHECRs. Any model that aims to the explanation of the problem of UHECRs 
must address all of the issues listed above (and possibly others) . 

In this section we consider in some more detail the issue of the GZK cutoff 
and the seriousness of its absence in the observed data. 

It is often believed that the identification of one or a class of nearby 
UHECR sources would explain the observations and in particular the absence 
of the GZK cutoff. This is not necessarily true. The (inverse of the) lifetime 
of a proton with energy E is plotted in fig. 1 (left panel) together with the 
derivative with respect to energy of the rate of energy losses b(E) (right panel) 
[the figure has been taken from ref. EDI)]. The flux per unit solid angle at 
energy E in some direction is proportional to tiq\{E)<&(E), where no is the 
density of sources (assumed constant), X(E) = c/((l/E)dE/dt) and $(J5) is 
the source spectrum. This rough estimate suggests that the ratio of detected 
fluxes (multiplied as usual by E 3 ), at energies E\ and E 2 is 



where in the last term we assumed that the source spectrum is a power law 
$(J3) ~ E-t. If for instance one takes E x = 10 19 eV (below E GZK ) and 
£ 2 = 3 x 10 20 eV (above E GZK ), from fig. 1 one obtains that K ~ 80 for 7 = 3 
and 7Z ~ 10 for 7 = 2.4. The ratio 1Z gives a rough estimate of the suppression 
factor at the GZK cutoff and its dependence on the spectrum of the source. 
For flat spectra (7 < 2) the cutoff is less significant, but it is more difficult to 
fit the low energy data H!3) (at E ~ 10 19 eV). Steeper spectra make the GZK 
cutoff more evident, although they allow an easier fit of the low energy data. 
The simple argument illustrated above can also be interpreted in an alternative 
way: if there is a local overdensity of sources by a factor ~ 1Z, the GZK cutoff is 
attenuated with respect to the case of homogeneous distribution of the sources. 
The question of whether we are located in such a large overdensity of sources 
was recently addressed, together with the propagation of UHECRs, in 
Assuming that the density of the (unknown) sources follows the density of 
galaxies in large scale structure surveys like PSCz HI) and Cfa2 HI), the 
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Figure 1: From Left panel) (\/E)dE/dt for a proton in IZ1) (curve 1), 

in HJ) (curve 2) and in HI) (black squares). The curve 3 is the contribution 
of the red shift. Right panel) The derivative db(E) / dE , with b(E) = dE/dt at 
z = 0. 



authors estimate the local overdensity on scales of several Mpc to be of order 
~ 2, too small to compensate for the energy losses of particles with energy 
above the threshold for photopion production. 

There is however another issue that the calculations in address, which 
is related to the statistical fluctuations induced by the process of photopion 
production. The large inelasticity of this process can be taken into account 
properly only through the use of Montecarlo calculations. When the Mon- 
tecarlo is applied to simulate the small statistics of events typical of current 
experiments, the fluctuations in the simulated fluxes above ~ 10 20 eV are very 
large, so that for flat spectra (~ E~ 2 ) the discrepancy between observations 
and simulations on the total number of events above ~ 10 20 eV is at the level 
of ~ 2cr (in agreement with the conclusions of Ref. HI)). The situation is 
represented in fig. 2 13) where the hatched regions show the uncertainties in 
the simulated fluxes. The data points are from AG AS A @> H). The bottom 
line of this section can be summarized in the following few points: 




Figure 2: Simulated fluxes with the AGASA statistics and"/ = 2.1. The sources 
are homogeneously distributed up to a maximum redshift z max = 1 . The differ- 
ent hatches refer to different cases of evolution of the sources H^) . 



1) the GZK cutoff is not avoided by finding sources of UHECRs that lie 
within the pathlength of photopion production, unless these sources are located 
only or predominantly nearby and are less abundant at large distances. 

2) The significance of the GZK feature depends on the fluctuations in the 
photopion production, and can be addressed properly only with a enhanced 
statistics of events with energy > 10 20 eV. 

4 The UHECRs engines 

Models for the origin of UHECRs can be strongly constrained on the basis on 
the criteria illustrated in the previous section. The challenge to conventional 
acceleration models, that are supposed to work at lower energy scales, induced 
an increasing interest for more exotic generation mechanisms, eventually re- 
quiring new particle physics. In this section the main ideas on production 
scenarios and their signatures will be summarized. 
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Figure 3: The Hillas Plot. 



4.1 Acceleration scenarios: the physics ol Zevatrons 

The accelerators able to reach maximum energies of order 1 ZeV have been 
named Zevatrons IH). The challenge for Zevatrons was recently discussed in 
detail by several authors E§1 EH): the main concept in this class of models is 
that the energy flux embedded in a macroscopic motion or in magnetic fields is 
partly converted into energy of a few very high energy particles. This is what 
happens for instance in shock acceleration. 

A discussion of all the models in the literature is not the purpose of 
the present paper, and in a sense we think it may not be very interesting. 
Nevertheless, it is instructive to understand at least which classes of models 
may have a chance to explain the acceleration to ZeV energies. In this respect, 
a pictorial way of proceeding is based on what is known as the Hillas plot H) . 
Our version of it is reported in fig. 4. 

In the Hillas plot the maximum energy is taken to be in its simplest form, 
as determined by the local magnetic field B and the size L of the accelerator: 
Emax = ZeBL. Here Ze is the electric charge of the accelerated particles. 
From fig. 4 it is evident that only 4 classes of sources have the potential to 
accelerate protons to ultra-high energies: 1) Neutron stars; 2) Radio Lobes; 



3) Active Galactic Nuclei; 4) Clusters of Galaxies. In the case of Iron, the 
situation becomes more promising for other sources, like the galactic halo or 
extreme white dwarfs. These sources would however have other problems that 
make them unlikely sources of UHECRs. 

The hillas plot does not include the effect of energy losses in the accel- 
eration sites. Photopion production limits the maximum energy achievable in 
clusters of galaxies to < a few 10 19 eV. These sources will therefore not be con- 
sidered any longer as sources of cosmic rays above the GZK energy. We also 
do not discuss here the so-called bursting sources. The prototypical example 
of these sources are gamma ray bursts, that have been proposed as sources of 
UHECRs EE]) . We refer the reader to recent literature treating this topic E3) . 

In the following we briefly summarize the situation with the other three 
classes of objects listed above. 

Neutron Stars 

The possibility that neutron stars may be accelerators of UHECRs was 
discussed in detail in Ref. HI) (and references therein). The main problem 
encountered in reaching the highest energies is related to the severe energy 
losses experienced by the particles in the acceleration site H' HI). Most of 
the mechanisms discussed in the literature refer to acceleration processes in 
the magnetosphere of the neutron star, where curvature radiation limits the 
maximum energy to a value much smaller that 10 20 eV. 

An alternative approach is to think of acceleration processes that occur 
outside the light cylinder of young neutron stars HI) . 

Rapidly rotating, newly formed neutron stars can induce the acceleration 
of iron nuclei through MHD winds outside the light cylinder HI). Although 
the mechanism through which the rotation energy of the star is converted into 
kinetic energy of the wind is not yet completely understood, it seems from 
the observations of the Crab nebula that a relativistic wind does indeed exist, 
with a Lorentz factor of ~ 10 7 EE'. Possible nuclei with charge Z26 = Z/26 
can be accelerated in young neutron stars to a maximum energy E max — 8 x 
10 20 Z 2 6-Bi3^3 fe , as estimated in Ref. 0). Here B13 is the surface magnetic 
field in units of 10 13 G and f^fc = il/3000s _1 is the rotation frequency of the 
star. Energies gradually smaller are produced while the star is spinning down, 
so that a spectrum ~ E^ 1 is produced by a neutron star. The process of escape 
of the accelerated particles becomes efficient about a year after the neutron star 



birth. Particles that are generated earlier cannot escape, but can produce high 
energy neutrinos in collisions with the ambient particles and photons 0) . The 
issue of the anisotropy due to the galactic disk is currently under investigation 

ID 

Active Galactic Nuclei 

Active galaxies are thought to be powered by the accretion of gas onto 
supermassive black holes. Acceleration of particles can occur in standing shocks 
in the infalling gas or by unipolar induction in the rotating magnetized accretion 
disk HI) . In the former scenario energy losses and size of the acceleration region 
are likely to limit the maximum energy of the accelerated particles to <C 10 20 
eV. In the latter case, the main limiting factor in reaching the highest energies is 
represented by curvature energy losses, that are particularly severe 0) unless 
moderately high magnetic fields can be kept with a small accretion rate. This 
may be the case of dormant supermassive black holes, possibly related to the 
so-called dark massive objects (DMO). Some 32 of these objects have been 
identified in a recent survey and 14 of them have been estimated to have 
the right features for acceleration of UHECRs . Had this model to be right, 
it would not be suprising that bright counterparts to the UHECR events were 
not found, since DMOs are in a quiescent stage of their evolution. 

Very little is known of DMOs as cosmic ray accelerators: the spectrum 
is not known, and neither is known their spatial distribution. It is therefore 
hard to say at present whether DMOs can satisfy the criteria listed in section 
|3|. In EH an interesting prediction was proposed: if UHECRs are accelerated 
by unipolar induction, they have to radiate part of their energy by synchrotron 
emission, resulting in the sources to become observable at TeV energies. 

Jets and lobes 

One of the most powerful sites for the acceleration of UHECRs is the 
termination shock of gigantic lobes in radio galaxies. Of particular interest are 
a subclass of these objects known as Fanaroff- Riley class II objects (FR-II), 
that can in principle accelerate protons to ~ 10 20 — 10 21 eV and explain the 
spectrum of UHECRs up to the GZK cutoff © . These objects are on average 
on cosmological distances. The accidental presence of a nearby source of FR-II 
type might explain the spectral shape above the GZK energy, but it would not 
be compatible with the observed anisotropy @' OH) . Nevertheless, it has been 
recently proposed that a nearby source in the Virgo cluster (for instance M87) 



and a suitable configuration of a magnetized wind around our own Galaxy 
might explain the spectrum and anisotropy at energies above ~ 10 20 eV 
as measured by AGASA. This conclusion depends quite sensibly on the choice 
of the geometry of the magnetic held in the wind. Several additional tests to 
confirm or disprove this model need to be carried out. 

4.2 The Top-Down Approach 

An alternative to acceleration scenarios is to generate UHECRs by the decay 
of very massive particles. In these particle physics inspired models the problem 
of reaching the maximum energies is solved by construction. The spectra of the 
particles generated in the decay are typically flatter than the astrophysical ones 
and their composition at the production point is dominated by gamma rays, 
although propagation effects can change the ratio of gamma rays to protons. 
The gamma rays generated at distances larger than the absorption length pro- 
duce a cascade at low energies (MeV-GeV) which represents a powerful tool 
to contrain TD models There are basically two ways of generating the 

very massive particles and make them decay at the present time: 1) trapping 
them inside topological defects; 2) making them quasi-stable (lifetime larger 
than the present age of the universe) in the early universe. We discuss these 
two possibilities separately in the next two sections. 

4.2.1 Topological Defects 

Symmetry breakings at particle physics level are responsible for the formation 
of cosmic topological defects (for a review see Ref . E§) ) . Topological defects as 
sources of UHECRs were first proposed in the pioneering work of Hill, Schramm 
and Walker E^. The general idea is that the stability of the defect can be 
locally broken by different types of processes (see below): this results in the 
false vacuum, trapped within the defect, to fall into the true vacuum, so that 
the gauge bosons of the field trapped in the defect acquire a mass mx and 
decay. 

Several topological defects have been studied in the literature: ordinary 
strings 0) , superconducting strings IH) , bound states of ma gne tic monopoles 



|51|i E3> , networks of monopoles and strings fc^D , necklaces u3 and vortons 
IH) . Only strings and necklaces will be considered here, while a more extended 
discussion can be found in more detailed reviews 0' H) . 



Ordinary strings 

Strings can generate UHECRs with energy less than mx ~ rj (the scale of 
symmetry breaking) if there are configurations in which microscopic or macro- 
scopic portions of strings annihilate. It was shown H) that self-intersection 
events provide a flux of UHECRs which is much smaller than required. The 
same conclusion holds for intercommutation between strings. 

The efficiency of the process can be enhanced by multiple loop fragmen- 
tation: as a nonintersecting closed loop oscillates and radiates its energy away, 
the loop configuration gradually changes. After the loop has lost a substantial 
part of its energy, it becomes likely to self-intersect and fragment into smaller 
and smaller loops, until the typical size of a loop becomes comparable with the 
string width 77 and the energy is radiated in the form of X-particles. Although 
the process of loop fragmentation is not well understood, some analytical ap- 
proximations HI) show that appreciable UHECR fluxes imply utterly large 
gamma ray cascade fluxes (see however 0) ) . 

Another way of liberating X-particles is through cusp annihilation H) , 
but the corresponding UHECR flux is far too low compared with 

observations. 

The idea that long strings lose energy mainly throu gh f ormation of closed 
loops was recently challenged in the simulations of Ref. ElP , which show that 
the string can produce X-particles directly and that this process dominates 
over the generation of closed loops. This new picture was recently questioned 
in Rcfs. H H. 

Even if the results of Ref. 0) are correct however, they cannot solve the 
problem of UHECRs 1^ : in fact the typical separation between two segments 
of a long string is comparable with the Hubble scale, so that UHECRs would 
be completely absorbed. If by accident a string is close to us (within a few 
tens Mpc) then the UHECR events would appear to come from a filamentary 
region of space, implying a large anisotropy which is not observed. Even if the 
UHECR particles do not reach us, the gamma ray cascade due to absorption of 
UHE gamma rays produced at large distances imposes limits on the efficiency 
of direct production of X-particles by strings. 

Necklaces 

Necklaces are formed when the following symmetry breaking pattern is 
realized: G — » H x U(\) — » H x Z^- In this case each monopole gets attached 
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Figure 4: Fluxes of UHECRs from necklaces HI) . 

to two strings (necklace) El) . The critical parameter that defines the dynamics 
of this network is the ratio r = mj [id where m is the monopole mass and 
d is the typical separation between monopoles (e.g. the length of a string 
segment). If the system evolves toward a state where r ^> 1, the distance 
between the monopoles decreases and in the end the monopoles annihilate, 
with the production of X-particles and their decay to UHECRs. The rate of 
generation of X-particles is easily found to be fix ~ r 2 /j,/t 3 mx- The quantity 
r 2 /i is upper limited by the cascade radiation, given by uj cas — \f-nf 2 [L = 
z/ttT 2 /i/ig (/ w i~ 0.5 — 1). The typical distance from the Earth at which 
the monopole-antimonopolc annihilations occur is comparable with the typical 

separation between necklaces, D ~ ( Jf J > 10(/i/10 6 Gel/ 2 ) 1/4 kpc. 

Clearly, necklaces provide an example in which the typical separation be- 
tween defects is smaller than the pathlength of gamma rays and protons at 
ultra-high energies. Hence necklaces behave like a homogeneous distribution of 
sources, so that the proton component has the usual GZK cutoff. This com- 
ponent dominates the UHECR flux up to ~ 10 20 eV, while at higher energies 
gamma rays take over. The fluxes obtained in Ref. HI) are reported in Fig. 
4, where the SUSY-QCD fragmentation functions lit were used. The dashed 



lines are for mx = 10 14 GeV, the dotted lines for mx = 10 15 GeV and the 
solid lines for mx = 10 16 GeV. The two curves for gamma rays refer to two 
different assumptions about the radio background at low frequencies 1^1) . 

4.2.2 Cosmological relic particles 

Super heavy particles with very long lifetime can be produced in the early 
universe and generate UHECRs at present U3 U3 In order to 

keep the same symbolism used in previous sections, we will call these particles 
X-particles. 

The simplest mechanism of production of X-particles in the early uni- 
verse is the gravitational production EH); particles are produced naturally in 
a time variable gravitational field or indeed in a generic time variable classical 
field. In the gravitational case no additional coupling is required (all particles 
interact gravitationally) . If the time variable field is the inflaton field <f>, a di- 
rect coupling of the X-particles to <f> is needed. In the gravitational production 
inflation is not required a priori, and indeed it reduces the effect. It can be 
shown that at time t, gravitational production can only generate X-particles 
with mass mx < H(t) < m^, where H(t) is the Hubble constant and is the 
inflaton mass. The authors in Ref. and 0) demostrated that the fraction 
of the critical mass contributed by X-particles with mx ~ 10 13 GeV produced 
gravitationally may be fix ~ 1. 

If the X-particles are directly coupled to the inflaton field, they can be 
effectively generated during preheating EH EH) . Alternative mechanisms for 
the production of X-particles are based on non-equilibrium thermal generation 
during the preheating stage 1^ . 

As mentioned in the beginning of this section, in order for X-particles to 
be useful dark matter candidates and generate UHECRs they need to be long 
lived (for a possible annihilation scenario see Ref. III). The gravitational 
coupling by itself induces a lifetime much shorter than the age of the universe 
for the range of masses which we are interested in. Therefore, in order to 
have long lifetimes, additional symmetries must be postulated: for instance 
discrete gauge symmetries can protect X-particles from decay, while being very 
weakly broken, perhaps by instanton effects HI) . These effects can allow decay 
times larger than the age of the universe, as shown in EH*. The slow decay 
of X-particles produces UHECRs. The interesting feature of this model is 



that X-particles cluster in the galactic halo, as cold dark matter ESP. Hence 
UHECRs are expected to be produced locally, with no absorption, and as a 
consequence the observed spectra are nearly identical to the emission spectra, 
and therefore dominated by gamma rays. The very flat spectra and the gamma 
ray composition are two of the signatures. Th e calculations of the expected 

The strongest signature of 
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fluxes have been performed in 
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the model is a slight anisotropy due to the asymmetric position of the sun in 



the Galaxy 
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V . More recently a detailed evaluation of the amplitude 
and phase of the first harmonic has been carried out in llJ and H). The 
two papers agree that the present data is still consistent with the anisotropy 
expected in the model of X-particles in the halo. Small scale anisotropics do 
not find an easy explanation in TD models, with possibly the exception of the 
SH relics H). 



5 Hints of New Physics? 

The possibility that at sufficiently high energies some deviations from known 
Physics may occur is of particular interest for UHECRs. Some attempts have 
been made to explain the events above 10 20 eV as a manifestation of some 
kind of Physics beyond the Standard Model of particle interactions. These 
suggestions became even more interesting after the recent claims for correlations 
of the arrival directions of UHECRs with ob ject s at large redshift. Some of 
these correlations are still subject of debate t3 . More recent studies result 
in a quite intriguing correlation with BL Lacs HI). Two of the BL Lacs 
in the sample used in HI) are in the error box of the two triplets of events 
detected by AGASA, and correspond to a distance of ~ 600 Mpc, much larger 
than the pathlength for photopi on p roduction of protons (however more than 
half of the objects used in Ref. E3^ have unknown redshift [Tkachev, private 
communication] ) . 

The absence of bright nearby counterparts to the UHECR events has first 
inspired theoretical proposals of neutrinos as primary particles, since these par- 
ticles are not affected by the presence of the CMB and can therefore propagate 
on cosmological distances without energy losses other than those due to the 
cosmic expansion. However, the small cross section of neutrinos makes them 
unlikely primaries: they simply fail to generate the observed showers in the 
atmosphere. There is one caveat in this argument: the neutrino nucleon cross 



section at center of mass energies above the electroweak (EW) scale has not 
been measured, so that the argument above is based on the extrapolation of 
the known cross sections, and simply limited by the weak unitarity bound H) . 
It has been proposed HOP that an increase in the number of degrees of freedom 
above the EW scale would imply the increase of the neutrino nucleon cross 
section above the standard model prediction. In particular, in the theories 
that predict unification of forces at ~ 1 TcV scale with large extra dimensions, 
introduced to solve the hierarchy problem 1^ , these additional degrees of free- 
dom arise naturally, and imply a neutrino-nucleon cross section that increases 
linearly with energy, reaching hadronic levels for neutrino energies of ~ 10 20 
eV (for string scale at 1 TeV). The calculations in Ref. 0^ show that the 
cross section remains too small to explain vertical showers in the atmosphere 
(for possible upper limits on this cross section see Ref. 

An alternative, even more radical proposal to avoid the prediction of a 
GZK cutoff in the flux of UHECRs consists in postulating a tiny violation of 
Lorentz invariance (LI). The main effect of this violation is that some processes 
may become kinematically forbidden 0). In particular, photon-photon pair 
production and photopion production may be affected by LI violation. The 
absence of the GZK cutoff would then result from the fact that the threshold 
for photopion production disappears and the process becomes kinematically 
not allowed. It is suggestive that the possibility of a small violation of LI has 
also been proposed to explain the apparent absence of an absorption cutoff in 
the TeV gamma ray emission from Markarian-like objects El) (see El) for a 
critical view of this possibility). 

6 Conclusions 

The increasing evidence for a flux of UHECRs exceeding the theoretical expec- 
tations for extragalactic sources has fueled interest in several models. These 
models aim to satisfy all the requirements imposed by observations on fluxes, 
spatial anisotropies and composition. At present, however there is no obvi- 
ous successful model. The situation might change with the availability, soon to 
come, of quite larger statistics of events, that will be achievable by e xpe riments 
like the Pierre Auger Project H) and EUSO/Airwatch/OWL H H). 

These future experimental efforts will be crucial mainly in three respects: 
1) the increase of statistics by a factor 100 for Auger and even more for the 



space based experiments will allow to strongly constrain theoretical models, 
and check whether the present excess is a (2 — 3)cr fluctuation or a physical 
effect. Moreover the small scale anisotropics, if real, will become stronger and a 
correlation function approach will definitely become appropriate to the analysis 
of the events; 2) the full sky coverage will finally allow a test of models based 
on local extragalactic sources and on galactic sources of UHECRs, through 
the measurement of the large scale anisotropy, that is currently spoiled by the 
limited spatial exposure; 3) a better determination of the composition of the 
UHECR events will represent a smoking gun either in favour of or against 
models: TD models would be ruled out if no gamma rays are found or if heavy 
nuclei represent the main component. On the other hand, iron-dominated 
composition would point toward a possible galactic origin, possibly related to 
neutron stars. 

The confirmation, on statistical grounds, of the association of UHECRs 
to distant cosmological objects like BL Lacs would represent a very strong 
indication of physics beyond the standard model. Either new interactions or 
some modification of fundamental physics would be needed to explain such a 
result. 
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